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1.  INT^< 


U  r-r  S  <  : 


ODUCTION 


The  search  for  a  Grand  Unified  Theory  in  quantum  physics  has  led  to  the  postulation  of  a  so-called 
(“fifth  foree^and  the  so-far-undetected  particle  associated  wth  propagation  of  this  field^^Concurrently, 
recent  geophysical  evidence  suggests  slight  deviation  of  measured  ^gravitatlonal^  attractipn  from 
Newton’s  inverse  square  formulation . ^^he  possibility  of  detecting  this  anomalous  “force’’  has  spawned 
the  ^ttjwer  experiment^ in  progress  by  the  dftr-Force  Geophysics  laboratory  ^AFGI^  In  this  experiment, 
deviation  of  measurements  made  at  various  heights  on  a  600-m  tower  from  values  predicted  by 
Newtonian  upward  continuation  of  ground-based  measurements  is  the  basis  for  interpretation. 


2.  OBJECTIVES  OF  THE  RESEARCH  EFFORT 

Though  departure  from  Newton's  Law  could  be  ascribed  to  various  phenomena  (for  example, 
spatial  variation  of  the  gravitational  constant),  the  field  is  modeled  assuming  an  anomalous  Yukawa 
contribution  consistent  with  a  potential  due  to  an  isolated  mass  of  the  form 


(1  +  ae-|ir) 
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(Chan  and  Paik6).  Objectives  are  to  analyze  field  behavior  to  determine  necessary  precision  of  measure¬ 
ment  and  to  study  effects  of  depth-related  and  lateral  inhomogeneities  using  models  simple  enough  to 
yield  such  information  unambiguously.  This  effort  is  designed  to  complement  simultaneous  efforts  in 
more  detailed  design,  analysis,  and  interpretation  of  the  "tower  experiment." 

3.  THE  FIELD  DUE  TO  A  UNIFORM  SPHERE 

Assume  a  sphere  of  uniform  density  p,  volume  V,  and  radius  R  centered  at  the  origin  (Figure  1).  Let 
p  ?  R  represent  the  distance  from  the  origin  to  a  point  P.  Let  us  find  the  Newtonian  and  Yukawa 
components  of  the  total  field,  denoted  by  Fn(P),  Fy(P),  and  F(P),  respectively. 

3-1  General  Approach 

For  each  element  dV  of  the  sphere,  let  lr  be  a  unit  vector  in  the  direction  from  P  to  dv,  let  r  be 
the  distance  from  P  to  dV,  and  let  lr  be  a  unit  vector  in  the  direction  from  P  to  the  origin.  Then  the 
field  can  be  written  as 

F(P)  =  FN(P)  +  FY(P) 

-  [Fn(P)  +  FY(P)]  1R 


where 


Fn(P)  =  P  Go 


in)  dv, 


(1) 


Fy(P)  = 


-a  p 


Or  i*)  dv- 


(2) 


i 

and  G0  lr  and  -a  G0  lr  are  the  gradients  of  the  Newtonian  and  Yukawa  potentials,  respec¬ 
tively,  due  to  each  element  dv. 

Specifically,  because  these  potentials  are  given  by 


4>n(P)  = 


G0pdV 

r 


(3) 


<MP) 


r 


(•*) 


6.  Chan,  HA.  and  Paik,  H.J.  (1984)  Experimental  Test  of  a  Spatial  Variation  of  the  Newtonian  Gravitational  Constant 
at  Large  Distances  in  Precision  Measurement  and  Fundamental  Constants  //,  Taylor,  B.N.  and  Phillips,  W.D.  (eds  ), 
Nad.  Bur.  Stand.  (U.S.),  Spec.  Publ.  617,  pp.  601 -606. 
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p 


Figure  1.  Notation  for  a  Uniform  Sphere 


their  gradients  (noting  that  we  have  chosen  lr  to  point  toward  dV)  are: 


grad  <}»N(P)  ~  °£  lr 

(5) 

grad  4>y(P)  =  -  e'^  (1  +  \ir)  ir. 

(6) 

Hence,  the  functions  f  and  g  are  given  by 

pH 

ll 

§ 

(7) 

g(r)  =  e"  (1  +  fir). 

a 

(8) 

In  what  follows,  let  h(r)  represent  f(r)  or  g(r),  as  appropriate. 

Our  problem  now  is  one  of  evaluating 

I(P)  =  /  1?  (ir  '  lfd  dV  (9) 

To  take  advantage  of  symmetry,  for  each  r  c  [p  -  R,  p  +  R]  let  S  represent  the  intersection  of  V  with  a 
spherical  shell,  of  radius  r  and  thickness  dr,  centered  on  P.  Let  (J>  =  cos-1  (lr  •  Ir).  and  for  each  d>, 
form  a  ring  in  S  of  radius  r  sin  4>  and  thickness  d<J>.  The  volume  of  the  ring  is  then  2-rr  (r  sin  4>)(rd4>)dr, 
the  projection  lr  Ir  for  each  dV  in  the  ring  is  cos  6,  and  r  is  constant  on  S;  hence 


3 


/P-R  /•*„ 

h(r)  /  sin  <(>  cos  4>  d4>  dr 

•I  Jo 


where  4>0  *»  sup  <j>.  Evaluating  the  inner  integral  yields 


/•r+R 

I(P)  *  tt  /  h(r)  sin2  <J>0  dr. 

»P-R 


(10) 


Using  the  law  of  cosines,  one  obtains 


sin2  <|>0 


.  / r2  +  p2  -  R2\2 

1  '  l — TZr— )  :  p  >  R’ 


(ii) 


1  -  (r/(2p))2 


;  p  =  R- 


3.2  The  Newtonian  Field 

As  a  check  on  the  above  method,  the  author  evaluated  the  above  integral,  using  h(r)  =  f(r)  =  1 . 
The  result  was 


I„(P)  -  *0-.  or 


(12) 


FN(P)  -  G0  p  i^4£, 
3  P 


(13) 


a  well-known  result  more  easily  obtained  via  the  Divergence  Theorem.  Note  that  in  the  case  of  a  uniform 
sphere,  the  Newtonian  field  outside  the  sphere  is  as  if  the  entire  mass  were  concentrated  at  the  center. 

3.3  The  Yukawa  Field 

For  the  Yukawa  field,  the  I(P)  integral  is  given  by 


/•P  +  R 

Iy(P)  =  ir  /  e_,ir  (1  +  pr)  sin2  4>0  dr. 

•'P-R 


(14) 


Using  Eq.  (11)  with  A2  =  p2  -  R2  gives 


IY(P)  -  tt  fj  *  e-^  (1  +  pr)  (1  -  A4/(4p2r2)  -  A2/(2p2)  -  r2/(4p2)J  dr. 
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The  Yukawa  field  does  not  have  a  closed-form  representation,  as  multiplying  out  the  integrand  gives  two 
terms  e "  “Vr  and  e_,ir/r2.  Slow  convergence  of  the  series  representations  rendered  this  approach 
computationally  infeasible  for  p  near  R. 

However,  the  integrand  of  Eq.  (14)  is  well-behaved,  consisting  of  the  product  of  a  function  dependent 
only  on  the  nature  of  the  field  and  the  purely  geometric  factor,  sin2  60  Rewriting  Eq.  (14)  in  terms  of  the 
dimensionless  parameter  r/S,  where  5  =  jjl  ~ 1  is  the  distance  from  an  isolated  point  mass  at  which  the 
Yukawa  potential  would  be  1/c  times  its  value  at  the  point,  yields 


Iy(P) 


= irs  r 


e-(rs>  (j  +  r/$)  gin2  <j>0  d(r/8). 


(15) 


where  sin2  <J>0  is  given  in  terms  of  r/5  by 


sin2  <t>0  = 


,  _  ((r/6)2  +  ((p-R)/6)((p  +  R).-8)\2 
'  ((p-R).fi  +  (p  +  R)/$j  Co'S)  '  ,P 


(p-R)8  +  (p+R)/8 


(16) 


Figures  2  and  3  show  plots  of  the  integrand  functions;  Figure  4  displays  IY  (r/5)/-n-8  as  a  function  of 
(p  -  R)/5.  An  interpretation  of  Figure  4  is  that,  for  a  given  8,  the  Yukawa  field  persists  to  larger 


Figure  2.  Plot  of  c  '  (1  +  v)  Versus  v 


r/R 


Figure  3  Sin  <i>0  Versus  r/R  for  Various  5/R  (Uniform  Sphere) 


distances  from  the  surface  as  the  radius  of  the  sphere  increases,  approaching  the  limiting  condition  of  a 
plane  surface  approximated  by  the  graph  corresponding  to  5/R  =  0  001. 


3.4  Yukawa  and  Newtonian  Fields  Compared 

According  to  the  definition  of  the  dimensionless  parameter  r/5  =  pr  for  the  Yukawa  field.  F.qs  (1) 
and  (2)  take  the  forms: 


Fn(P) 


p  Go  Tt  S 


(1>  sin2  <b0  d(r'5). 


(17) 
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Jy(P)  -  Fv(P)/(ap  G0  tt  5), 


(20) 


corresponding  to  the  values  of  the  respective  integrals.  Eqs.  (17)  and  (13)  together  give 


Jn(P)  - 


4 

3(p/R)2  (8/R) 


(21) 


corresponding  to  the  Newtonian  field. 

Perhaps  the  most  obvious  comparison  of  the  Yukawa  and  Newtonian  field  components  is  via  their 
ratio,  specifically,  by 


1  Fv(P)  _  M2- 

«  Fn(P)  Jn(P)  ’ 


as  depicted  in  Figure  5.  The  range  of  the  Yukawa  field  compared  to  the  Newtonian  field  is  limited;  the 
ratio  Jy/Jn  for  each  curve  drops  by  three  decades  from  its  surface  value  within  10  8  units  of  distance 
from  the  surface.  Rollover  for  each  curve  begins  at  or  before  (p  -  R)/8  =  1,  but  in  contrast  to  Figure  4, 
the  smaller  8/R  values  correspond  to  slightly  earlier  and  more  rounded  rollovers.  This  behavior  is 
explained  by  Jn  being  more  nearly  constant  near  the  surface  for  larger  values  of  R/S,  as  the  surface 
corresponds  to  points  fardier  out  on  the  inverse  square  dependence.  The  main  feature  of  Figure  5  is  the 
overall  decrease  in  the  height  of  the  curves,  in  particular  the  surface  values  of  Jy/Jn,  as  6/R  decreases. 
Clearly,  the  JN  becomes  much  greater  than  JY  as  R  increases,  because  JN  is  proportional  to  the  volume 
whereas  Jy  is  effectively  more  locally  defined  for  R  >  >  8. 

Because  Fy  and  FN  cannot  be  measured  individually  at  present,  interpretation  from  actual  measure¬ 
ments  must  be  based  on  some  other  parameter,  definable  in  terms  of  the  total  field.  Such  a  parameter 
will  be  presented  shortly. 

4.  THE  FIELD  DUE  TO  A  LAYERED  SPHERE 

Consider  a  set  of  concentric  spheres  Sj,  S2,...,  Sn  of  radii  R)  >  Rj  .  .  .  >  RN  >  0,  respectively. 
Define 


Pi  *  density  of 


shell  Sj  D  Sj+ 1;  i  =  1,  2 . N-l, 


sphere  Sn  ;  i  =  N. 


Because 


Pi 


Pi  ;  i  =  1 

(Pi  ~  P1-1)  +  Pi  - 1  -  >  =  2 . N, 


(23) 
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Figure  5.  Jy/Jn  Versus  (p  -  R)/5  for  Various  8/R  (Uniform  Sphere) 


the  density  distribution  can  as  well  be  duplicated  by  superposing  the  spheres  themselves,  haring 
corresponding  densities 


A,  = 


Pi  ;  >  =  1 

Pi  ~  Pi-i  ;  i  =  2 . N, 


(24) 


The  field  due  to  a  layered  sphere  is  then  the  superposition  of  fields  due  to  individual  spheres  whose 
radii  correspond  to  the  layer  boundaries  and  whose  densities  are  given  in  terms  of  the  layer  densities  by 


Eq.  (24).  Specifically,  if  Jn  and  Jy  are  defined  as  in  Eqs.  (19)  and  (20),  where  now  FN  and  FY  are  for  the 
layered  case  and  pi  replaces  p  in  the  definition,  one  obtains 


J  *  (P) 


i 


l-l 


YiJ  *i(P,Si) 


(25) 


where 


Yi  =  Ai/pi, 

*  *  N  or  Y, 


and  J  *j  (P,  Sf)  represents  the  value  of  J*  at  point  P  due  to  a  single  sphere  of  density  A;  According  to 
the  previous  section,  the  J*t  (P,  S|)  quantities  are  given  by 


JNi(P >  SO 


4 

3(P/R02  (8/RO’ 


Jyi(P,S  0 


r  +  R. 


(1  +  v)  sin2  4>0j  dv 


(26) 


(27) 


where 


sin2  4>oi 


v2  +  «p-R)/S]  f(p  +  R)/3]2 
v[(p-R)/S  +  (p  +■  R)/S]  ’  p 

_ V _ 2  n  =  R 

(p-R)/S  +  (p  +  Rj/'S  ’  P 


(28) 


5.  SENSING  THE  YUKAWA  COMPONENT 

This  section  defines  a  parameter  sensitive  to  the  Yukawa  field.  The  effects  of  changing  model  and 
field  parameters  are  then  treated.  Finally,  the  topic  of  required  measurement  accuracy  is  addressed. 
Layered  or  uniform  sphere  models  are  assumed  throughout  this  section. 


5.1  Interpretation  Parameters 

A  primary  reason  for  the  layered  sphere  model  is  to  enable  study  of  the  effects  of  inhomogeneities 
on  the  interpretation  parameters  available  from  field  measurements.  Hence,  let  us  define  such  parame¬ 
ters  according  to  the  criteria  that  the  parameter  must  be  (1)  calculable  from  field  data,  and  (2)  sensitive 
to  the  existence  of  the  Yukawa  field. 
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If  Go,  pi,  and  o  were  known  percisely,  a  parameter  such  as  F£4(P)  called  "apparent  Newtonian 
surface  field,”  defined  as 


F£,(P)  *  the  value  at  the  surface  of  a 
Newtonian  field  having  the  observed  value  F(P)  at 

point  P,  given  the  model,  (29) 


would  contain  all  the  information  available  from  the  measured  total  field.  If  this  parameter  is  "nor¬ 
malized”  by  dividing  by  the  total  surface  field  actually  observed,  this  new  parameter  has  the  advantage  of 
being  immune  to  inaccuracies  in  the  gravitational  constant  and  the  first  layer  density.  Mathematically,  this 
parameter  is  described  by 


FS,  „  F(P)  (p/R)2 
"W  F° 


(30) 


in  terms  of  observed  field  values,  and  by 


Un(P)  ~  «  Jy(P)]  (P/ft)2 

Lft  "  «J?1 


(31) 


in  terms  of  the  JN  and  Jy  integrals.  If  the  field  were  totally  Newtonian,  this  parameter  would  be  1  for  all 
p  &  R,  as  the  numerator  of  Eq.  (30)  then  is  exactly  *N.  according  to  Eq.  (21).  Hence,  deviation  of  this 
parameter  from  unity  is  indicative  of  non-Newtonian  behavior,  assuming  perfect  measurements  over  the 
Earth  of  the  model. 

Because  departure  from  “1”  is  of  interest,  a  new  parameter,  equivalent  to  this  departure,  is  what 
will  be  called  the  <i»  parameter: 

+(P)  -  T?  -  1  (32) 


In  terms  of  the  Jn  and  Jy  integrals,  one  obtains 


*l»(P)  = 


-«  Uy(P)(P/R)2  ~  J?) 
JS  -  «J? 


(33) 


by  substitution  from  Eq.  (30),  the  key  step  being  to  note  that  Jn(P)(P/R)2  ~  J\.  which  arises  in  the 
numerator,  is  zero. 

A  useful  feature  of  the  parameter  is  that  for  h/ R  <  <  1  (as  for  the  Earth)  and  for  a  <  0.1,  aj°  is 
much  less  than  in  the  denominator.  Thus,  in  practical  cases,  il»  can  be  studied  via  »b/a.  which  is 
virtually  a-independent  (see  Table  1). 
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Table  1.  iji/a  Versus  (p  -  R)/S  for  a  =  0.004, 
0.1  (8/R  =  1/10000) 


(P  -  R>/« 

S/a  (s  =  0.004) 

S/6  (a  =  0.1) 

0.0095 

1.42359  x  10~6 

1.43261  x  10“6 

0.0560 

8.17235  x  10"6 

8.17246  x  10"6 

0.3147 

4.04967  x  10" 5 

4.04973  x  10" 5 

1.7798 

1.24694  x  10" 4 

1.24696  x  10-3 

10.0624 

1.49994  x  10" } 

1.49996  x  10"} 

5.2  Influence  of  Inhomogeneity  With  Depth 

Consider  a  two-layer  sphere  of  radius  Rj.  Changes  in  the  first-layer  thickness  Hj  =  Rj  -  R2  and  in 
the  resistivity  contrast  p2/pj  clearly  change  the  scale,  but  not  the  shape  of  the  Newtonian  dependence 
for  p  >  Ri,  as  the  Newtonian  field  for  p  >  Rj  is  the  sum  of  two  terms  having  1  /p2  as  a  common  factor. 
Considering  the  fact  that  the  density  contrast  does  not  appear  inside  the  Jv  integral  of  Eq.  (27),  it 
follows  that,  for  fixed  8/Ri  and  (Rj  -  Rj)/8,  variations  in  P2/P1  affect  only  the  vertical  positioning  and 
not  the  shape  of  the  ’if/a.  curve  (Figure  6).  Comparing  these  curves  (for  6/Rj  =  1/10000)  with  those  of 
Figure  7  (for  8/Rj  =  1)  shows  that  this  dependence  on  p2/pi  tends  to  disappear  as  Rj  becomes  less 
than  a  few  multiples  of  8.  Also  note  from  Figure  6  that  the  p2/p i  dependence  tends  to  disappear  as  Rt 
becomes  less  than  a  few  multiples  of  8.  Also  note  from  Figure  6  that  the  p2/pi  dependence  of  curve 
position  decreases  as  p2/pi  becomes  greater  than  1. 

Less  mathematically  obvious,  but  computationally  verified  to  four  significant  figures  for  the  curves  of 
Figure  8,  is  that  for  fixed  8/Ri  and  p/2/pi,  the  shape  of  the  iji/a  curve  is  invariant  with  respect  to 
first-layer  thickness.  Again,  the  R2/Ri  dependence  is  subdued  for  R  less  than  a  few'  multiples  of  8.  Note 
from  Figure  8  that,  for  fixed  (P  -  Ri)/8  and  Rj  >  >  8,  log(t)//5)  varies  nearly  linearly  with  (Rj  - 
R2)/8  for  (R2  ~  Ri/8  <  1  whereas  for  (R2  -  Rp/j  >  4  iog(*/a’i)  is  virtually  independent  of  variations  in  this 
parameter. 

In  summary,  then,  layer  thickness  and  density  variations  with  depth  affect  only  the  relative  magni¬ 
tudes  of  the  Newtonian  and  Yukawa  components  on  or  outside  the  sphere.  Each  component  has  its  ow-n 
falloff-with-distance  characteristic,  whose  shape  is  invariant  with  respect  to  layer  parameters,  given  a 
specific  outer  radius  R2  and  the  field  parameters  a  and  8.  For  models  of  practical  interest  (8/R  <  <  1), 
the  above  is  essentially  true  independent  of  a,  which  effectively  enters  only  as  a  scaling  factor. 

As  a  specific  example,  Figure  9  shows  t| i/a  curves  for  two-layer  models,  where  Rj  is  the  radius  of  the 
Earth  (6371  km)  and  8  =  200  m.  The  curves  all  have  the  same  shape.  Thus,  if  the  Earth  were  a  layered 
sphere,  measuring  its  field  and  calculating  the  ifr  parameter  at  various  heights  would  ideally  give  either  a 
curve  similar  to  Figure  9  (scaled  by  a  factor  of  a),  indicating  presence  of  a  nonzero  Yukawa  component, 
or  else  essentially  zero  values  for  the  parameter,  indicating  that  any  Yukawa  component  present  was  too 
small  to  be  reliably  detected. 

The  face  that  the  curves  roll  over  around  H  =  600  m  means  that,  if  8  is  about  200  m,  measurements 
to  a  height  of  about  600  m  should  be  made  to  ensure  seeing  the  characteristic  features  of  the  curve 
(ramp  followed  by  rollover).  That  the  curves  for  various  values  of  hi  are  spread  only  over  a  range  of  a 
factor  of  2  suggests  that  ability  to  resolve  the  curve  is  affected  little  by  first-layer  thickness.  A  similar 
argument  applies  to  density. 

Note  that  the  effects  of  lateral  inhomogeneities  have  not  yet  been  considered.  Also,  curves  for 
multi-layer  spheres  have  the  same  shape  as  those  for  two-  and  one-layer  models. 
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Figure  6.  t|//a  Versus  (p  -  Ri)/B  for  Various  P2/P1  (Hi/8  *  1/4,  Rj/8  =  1/10000) 


5.3  Influence  of  Field  Parameters  a  and  S 

From  the  previous  discussion,  it  follows  that  the  shape  of  the  log-log  plot  of  4<  versus  h  is 
determined  completely  by  the  value  of  8/R.  For  8/R  <  <  1: 

1.  If  R  is  fixed,  lateral  positioning  of  the  curve  of  4*  versus  h  is  determined  solely  by  8.  Specifically,  the 
rollover  portion  of  the  curve  occurs  around  0.68  <  h  <  58. 

2.  That  4»/ot  is  virtually  independent  of  a  allows  the  effects  of  a  and  8  on  vertical  positioning  to  be 
studied  independently  of  each  other.  Clearly,  if  t|>/a  is  independent  of  a,  then  the  value  <l»a 
appproached  by  4/  as  h  increases  is  proportional  to  a.  If  a  is  held  fixed,  t)>a  turns  out  to  be 
proportional  to  8. 


For  example,  for  a  uniform  Earth  of  radius  6371  km,  and  for  8  set  to  1,  v|/a/ot  has  the  value 
2.36  X  10~*.  Hence,  for  this  Earth  and  any  a  and  8,  one  can  calculate 


<|ta  =  a  8  (2.36  x  lO"*). 


(34) 
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Figure  7.  ty/a  Versus  (p  -  Ri)/§  for  Various  P2/pi(Hj/8  =  1/4,  Ri/8  *  1) 


Because  the  influence  of  layer  parameters  on  vertical  positioning  of  the  ill  curve  is  relatively  minor,  the 
above  equation  will  be  a  useful  guideline  in  determining  required  accuracy  of  measurement  in  the  next 
section. 

A  clear  implication  of  (1)  and  (2)  above  is  that  8  could  theoretically  be  determined  (by  curve- 
matching  or  geometric  construction,  for  example)  from  horizontal  positioning  of  the  rollover.  Knowing  8, 
one  could  then  obtain  a  according  to 


=  a  8  »|j*, 


where  tjt*  is  a  reference  value  for  the  particular  model  type.  However,  it  seems  likely  that  the  estimate  of 
8  from  a  measurement-derived  curve  would  be  much  more  prone  to  error  than  the  a8  product  due  to 
uncertainty  in  curve  location  associated  with  measurement  errors.  Again,  geologic  noise  due  to  deviation 
of  the  real  Earth  from  that  of  the  model  has  not  yet  been  considered. 
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* 

Figure  8.  >|>/a  Versus  (p  —  R)/8  for  Various  Hi/8  (P2/PJ  =  4,  Rj/S  =  1/10000) 


5.4  Required  Accuracy  of  Measurement 

Assume  measurements  of  the  total  field  are  to  be  taken  at  various  altitudes  above  the  surface  of  a 
layered  sphere  with  8/R  <  <  1.  How  is  inaccuracy  in  the  measurements  related  to  resolvability  of  the 
^•parameter  curve  obtained? 

Assume  that  possible  inaccuracy  of  measurement  is  independent  of  height,  so  that  the  measured 
value  Ft  of  the  field  at  point  P  at  or  above  the  surface  is  described  by 


Ft  *  Ft  ±  e 


(35) 


Then,  according  to 


Ft  (p/R)2 

F? 


1 


the  calculated  value  of  4i  at  P  satisfies 


(36) 
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For  measurements  within  a  few  kilometers  of  the  surface  of  the  Earth,  p/R  »  1  and  Ft/F£  =  1:  thus, 
one  obtains 


+  (40) 

F?  F? 

Next,  consider  the  log-log  plot  of  as  a  function  of  height  h  (the  solid  curve  of  Figure  10).  According  to 
inequality  (40),  pairs  (hit  4»i)  calculated  from  measurements  will  lie  in  the  envelope  between  curves  1 
and-  2,  which  correspond  to  the  right  and  left  sides  of  the  inequality.  Note  that,  because  i|>  -*  0  as  h  -* 

0,  there  exists  some  minimum  value  of  h,  shown  as  hc,  at  which  the  value  of  >1)  becomes  smaller  than  the 
possible  error  2c/F£  in  its  estimated  value.  However,  points  to  the  right  of  hc  should  contain  more 
information  about  the  location  of  the  “actual”  curve.  How  well  the  curve  is  likely  to  be  approximated 
depends  on  the  tightness  of  the  envelope  between  curves  1  and  2  and  the  length  of  its  leftmost  extent, 
determined  by  2e/F£. 

Because  all  curves  are  the  same  shape,  all  have  rollovers  that  occur  at  the  same  multiple  of  5,  and 

the  asymptotic  value  r|»a  =  lim  4>  depends  primarily  on  the  parameters  a  and  5,  constraints  on 
n-»«j 

accuracy  can  be  formulated  largely  in  terms  of  plot  geometry. 

To  begin,  ramp  slopes  are  very  near  45  degrees,  and  all  curves  are  shaped  such  that  a  vertical  line 
h  =  58  lies  to  the  right  of  the  rollover  and  intersects  the  extrapolated  ramp  a  factor  of  5  higher.  A 
reasonable  criterion  for  acceptable  curve  resolution  would  be  to  specify  factors  u  and  v  such  that 


<|»/|>  -  (2e/F!J)l  <;  v  for  55/u  *  h  s  58;  (41) 

specifically, 

i|»/[i|i  -  (2e/F^)]  =  v  when  h  =  58/u,  (42) 

in  view  of  the  fact  that  the  envelope  between  curves  1  and  2  widens  toward  the  left. 

From  Figure  10,  one  obtains 


<|»  -  (2e/F$  pMi)  at  h  -  58/u; 


(43) 


5s 

2  u 


Figure  10.  i|/-Curve  Resolution  Determination  (Curve  1:  Ordinate  =  tli  +  2«/F?: 
Curve  2:  Ordinate  =  i|>  -  2e/F^) 


hence  the  constraint  (42)  becomes 

v  =  [58/ (uhc) ] /| 58/ uhc)  -  1] 


and  furthermore, 


5*./(2e/F?)  -  58/hc. 


(45) 


Eliminating  he  from  the  above  equations  and  solving  for  e  gives 


e 


5  (v  -  1) 
2  u  v 


Ft- 


(46) 


For  example,  suppose  R  =  6371  km,  and  S  and  a  are  postulated  as  0.2  km  and  0.008,  respectively. 
Assuming  a  uniform  Earth,  one  obtains  from  Eq.  (34) 


-  (0.008) (0.2) (2. 36  x  10“4)  =  3.76  x  10"'. 


Suppose  it  is  desired  that,  over  the  two  decades  of  h  to  the  left  of  h  =  58,  estimated  values  of  are  not 
to  depart  from  <1>  by  more  than  a  factor  of  1.2.  Then  u  =  100,  v  =  1.2,  and 

5(3-76  x  1Q-7)(1.2  -  1) 
v  r  2(100)(1.2) 

=■  1.57  x  IQ"9. 


Thus,  if  «*  1000  cm/s2,  the  tolerable  error  in  measurement  of  Ft  is 


e  »  1.6  x  10~6  cm/s2. 


(47) 


5.5  Atmospheric  Corrections 

Suppose  our  spherical  layered  Earth  has  an  atmosphere.  If  a  measurement  is  made  at  height  h 
above  the  surface,  what  corrections  should  be  applied  to  account  for  the  field  due  to  the  air? 

To  correct  the  Newtonian  component  to  the  value  that  would  have  been  measured  if  no  atmosphere 
were  present,  one  simply  subtracts  the  field  due  to  the  portion  of  the  atmosphere  that  lies  below  the 
alititude  of  the  measurement.  That  this  procedure  works  is  due  to  model  symmetry  and  to  the  fact  that, 
at  each  point  in  space,  the  divergence  of  Fn  is  proportional  to  the  mass  density  there.  Hence,  the 
volume  integral  in  the  Divergence  Theorem  is  essentially  the  mass  of  the  volume. 

However,  div  Fy  can  be  nonzero  even  in  free  space.  Thus,  if  points  at  which  div  Fy  =  0  are  called 
“sources"  of  the  Yukawa  field,  one  could  say  that  these  sources  are  "smeared”  spatially,  relative  to  the 
density  distributions  associated  with  them.  An  implication  is  that  the  atmosphere  both  above  and  below 
the  point  in  question  must  be  considered  in  calculating  the  Yukawa  field  due  to  the  atmosphere. 

Suppose  the  required  precision  in  measurement  (say  F  =  F  ±  «)  is  obtainable  and  the  Newtonian 
atmosphere  correction  is  made.  Whether  a  Yukawa  atmosphere  correction  is  necessary  depends  on 
whether 


Fy  air  ^  ^  Fy  Eartht  and  (48) 

Fy  air  <  <  «  (49) 


To  answer  the  first  criterion,  let  us  construct  a  crude  model  of  the  atmosphere  and  calculate 

Fy  air/Fy  Earth  at  a  height  of  0.6  km  (if  our  measurements  were  on  a  0.6-km  tower,  the  specified  quantity 

would  be  largest  at  the  top).  Using  p  =  kP/T  and  p0  =  0.00129  g/cm3  gives 


p  =  0.391  P/T 


(50) 


where  p  -  air  density  (pc  at  sea  level),  P  =  pressure  (atm.),  and  T  =  temperature  (deg  K).  Using 
dT/dh  =  -  6.5  deg  K/km  and  P  =  exp  [  -  h  (cm)/(2930  T)]  (Allen"’),  with  T  =  303  deg  K  gives 


p(h)  =  0.0013  -  1.18  x  10- 4  h  g/cm3, 


(51) 


with  h  in  kilometers,  as  a  linear  approximation  for  air  density  as  a  function  of  height.  Using  5  =  0.2  km 
and  layer  boundaries  at  0.6  ±  2n  (8/8)  (stopping  above  at  0.6  +  88  and  below  at  the  surface)  and 
assigning  to  each  layer  its  mean  air  density  gives 


Jy  =  1.92978,  Jy  =  1.95857 


corresponding  to  the  air  above  and  below  0.6  km,  respectively.  A  uniform  F.arth  has  Jy  Earth  =  0.098  at 
h  =  0.6  km.  Then,  using 


Fy  ajr  Pi  Jy  +  Pi  Jy 
Fy  can!  Pi  Jy  earth 


where  p  =  0.001231  g/cm3  and  pi  =  0.001228  g/cm3  are  the  densities  of  the  first  layers  above  and 
below  0.6  km,  respectively,  and  pi  =  2.57,  one  obtains 


Fy  air/Fy  earth  m  117  X  10  3 


(52) 


7.  Allen.  C.W.  (1981)  Astropbysical  Quantities.  3rd  ed.,  Univ.  London,  p.  120. 
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at  h  -  0.6  km.  Thus,  it  appears  that  no  correction  for  the  Yukawa  field  due  to  the  atmosphere  is 
necessitated  by  the  first  criterion. 

For  the  second  criterion,  Fy  air  <  <  e,  let  us  calculate  Fy  air  at  the  surface,  at  which  its  value  is 
largest.  Using  a  similar  layered  air  model  gives  J°  air  =  1.9618.  The  total  field  is  about  1000  cm/s2  and 
/n  E«th  *  4  274  *  lo4;  hence, 


Pi  air  Jyair 

F£  «  1000  a  - - -  =  1.8  x  10'7  cm/s2, 

Pearth  Jsjearth 


(53) 


where  F°  is  approximated  by  F£.  This  value  is  about  a  factor  of  10  less  than  e  =  1.6  x  10' 6  cm/s2  of 
the  previous  section,  so  for  these  examples,  one  need  not  correct  for  the  Yukawa  component  due  to  the 
atmosphere. 

6.  EFFECTS  OF  LATERAL  INHOMOGENEITIES 

The  extent  to  which  lateral  variations  in  density  affect  reliable  sensing  of  the  Yukawa  field  will  now 
be  considered.  To  enable  basic  understanding  of  field  behavior  in  terms  of  a  few  model  parameters,  the 
models  are  deliberately  simple. 


6.1  The  Field  Due  to  a  Uniform  Rod 

Assume  a  point  P  whose  distances  from  the  center  of  a  rod  are  q  units  measured  perpendicular  to 
the  rod  and  p  units  measured  parallel  to  the  rod  (Figure  11).  Let  the  rod  be  of  length  L  and  be 
described  by  its  linear  density  <r  =  pA,  where  A  is  its  cross-sectional  area.  The  Yukawa  component 
parallel  to  the  rod  is 


Fv,(P) 


-a  00  o’  / 


(54) 


where 


cos  4>  -  Vl  -  (q/r)2, 


(55) 


and  <1>  is  as  shown  in  Figure  11.  letting  v  =  pr  =  [rl)m/5  then  yields 


Kyj  (P) 


-  r 


vjp+iu)*/i 


(L+v) 


dv, 


V(p-L/2)z  +  q2/fc 


(56) 


where 
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Figure  11.  Yukawa  Components  of  Field  due  to  a  Rod 


Fy#(P)  -  Kyj  (P). 


The  component  perpendicular  to  the  rod  is  given  by 


„  „  A2  e~^r  (1  +  jir)  sin  <t>^, 

FYjl(P)  =  —a  G0  n  J  - p - dz> 


where  sin  4>  =  q/r.  The  substitution  v  =  r/S  yields 


KYl  (P)  =  (q/8)  f  ' 


(p-l •  e-v  (j  + 

v2  Vv2  -  (q/8)2 


\  (p-b2)-  +  q2  h 


where 


Fyx  (P)  =  ~ gG°--  KYi  (P). 


The  corresponding  Newtonian  components  are  described  by 


KNj  =  (B  -  A)/AB, 
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Kn. 


A  VB2  -  Q2  -  BVA2  -  Q1 
Q  A  B 


where  A  =  V(p  -  L/2)2  +  q^/8,  B  =  V(p  +  L/2)2  +  q2/S,  and  Q  =  q/5. 


(62) 


6.2  A  Rectangular  Block  In  a  Layered  Medium 

Consider  a  rectangular  parallelepiped  of  uniform  density  po  imbedded  in  a  much  larger  layered 
sphere  such  that  the  “top"  of  the  parallelepiped  is  essentially  parallel  to  the  locally  nearly  planar  surface 
(Figure  12). 

Subdivide  the  block  into  n  identical,  vertically  oriented  parallelepipeds.  Further  subdivide  each  smaller 
block  into  m  sections  defined  by  the  intersection  of  the  block  with  various  layers.  If  each  section  is 
assigned  an  appropriate  density  difference  and  then  replaced  by  a  vertical  rod  of  similar  height  and  mass 
passing  through  the  section's  center,  the  Yukawa  field  at  point  P  (vertical  component),  given  by 

Fy(P)  -  -  a  G0  p,  it  5  JJyl  (p)  +  Jyo(P)],  (63) 
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Figure  12.  >1>  Curves  in  Presence  of  Offset  Block 
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requires 


jv  m  -  s  i  h  (p>. 


i-l 


Pi 


i«  1 


(64) 


where  JyL  and  Jy0  are  for  the  field  of  the  layered  sphere  and  for  that  of  the  block  in  excess  of  that  of  the 
sphere,  respectively,  A  is  the  cross-sectional  area  of  each  section,  pj  is  the  density  of  the  layer  intersected  by 
thejth  section  of  a  block,  and  Kym  is  the  expression  of  Eq.  (56)  calculated  for  the  j,h  section  of  the  i'h  rod. 
The  expression  for  JNjj(P)  is  similar  except  that  KNij(P),  corresponding  to  Kq.  (5").  appears  in  place  of  KY||(P) 


6.3  The  Parameter  t|i 

The  parameter  t|i,  defined  previously  consistent  with 


Ft  (p/R)2 

F? 


-  1, 


is  calculated  for  the  block-and-sphcre  model  as 

Jn„(P/R)2  -  Jn0  -  a  Or  (p/R)2  ~  J?) 
*  ”  Jn  ~  «Jy 


(65) 


(66) 


where 


J.  =  J*L  +  J.,„ 

L  designates  due  to  the  layered  sphere, 

D  designates  due  to  density  difference  between  block  and  sphere, 
*  is  “N”  for  Newtonian  or  “Y"  for  Yukawa  component,  and 
superscript  “o”  indicates  surface  value. 


Lateral  inhomogeneity  distorts  the  location  and  shape  of  the  t|»  curve  from  that  of  a  layered  Earth.  The 
models  of  Figure  12  are  crudely  similar  to  the  case  of  a  broad  hill  near  the  measurement  site.  For  a  = 
0.008,  the  typical  rollover  is  masked  by  the  occurrence  of  higher  tli  values  near  the  surface,  relative  to  those 
for  the  homogeneous  model.  These  curves  lie  only  a  factor  of  1.5  above  those  for  a  =  0  (no  Yukawa  field) 
to  the  left  of  the  rollover;  hence,  for  small  a,  the  shape  of  the  curve  can  be  primarily  a  geometric  effect 
caused  by  distortion  of  the  Newtonian  field  near  the  inhomogeneity.  The  curve  for  a  =  0.08  shows 
behavior  more  typical  of  a  layered  model  and  lies  everywhere  at  least  a  factor  of  6  above  that  for  a  =  0 
due  to  greater  strength  of  the  Yukawa  component. 

These  results  indicate  that  a  reliable  interpretation  of  an  experiment  using  the  4>  parameter  requires 
that  the  portion  of  the  field  due  to  lateral  inhomogeneities  be  subtracted  before  the  <]>  curve  is  calculated. 
An  approach  that  avoids  using  a  parameter  defined  in  terms  of  one  specific  model  type  would  consist 
essentially  of  comparing  measurements  at  various  heights  with  the  upward-continued  surface  field, 
assumed  to  be  Newtonian. 
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7.  CONCLUSIONS 


For  a  sphere  or  layered  sphere  model,  the  shape  of  the  «]/-  parameter  curve  is  determined  by  the 
values  of  8/R  and  a;  changes  in  density  or  layer  thickness  affect  only  the  vertical  positioning  of  the 
curve.  For  8/R  <  <  1,  as  for  the  Earth,  the  shape  of  the  i|»  curve  is  virtually  a-independent  as  well,  as 
a  enters  into  essentially  as  a  scaling  factor.  Hence,  such  models  enable  study  of  the  relationship  between 
»|i  curve  discemibility  (Yukawa  field  detection)  and  measurement  precision,  exclusive  of  effects  of  lateral 
inhomogeneities. 

The  major  limitation  of  these  models  is  that  the  Earth  is  not  free  from  lateral  density  and  thickness 
variations.  Thus,  although  the  Yukawa  field  is  a  rather  local  effect,  distortion  of  the  much  larger  Newtonian 
component  due  to  lateral  density  variations  can  render  the  curve  uninterpretable  in 
raw  form. 

The  effect  of  atmospheric  mass  on  the  Yukawa  component  appears  negligible,  and  the  effect  on  the 
Newtonian  component  is  removed  easily. 


8.  RECOMMENDATIONS 

One  might  explore  the  degree  to  which  the  ijj-parameter  approach  in  the  presence  of  lateral 
inhomogeneities  could  be  salvaged.  In  the  experimental  environment,  such  an  approach  might  involve 
modeling  of  the  subsurface  followed  by  determination  of  the  residual  field  in  excess  of  that  of  a  "best-fit" 
spherically  layered  model.  Subtracting  the  residual  field  at  measurement  points  on  the  tower  should  then 
yield  a  vertical  field  profile  as  if  no  lateral  variations  were  present.  A  major  concern  with  this  method  is 
that  the  modeling  phase  (which  must  assume  a  and  8)  could  result  in  a  mass  distribution  that  would  yield 
erroneous  corrections  for  tower  measurements  while  still  reproducing  the  ground  measurements.  This 
phenomenon  would  be  due  to  the  Newtonian  and  Yukawa  fields  having  different  ranges  of  influence. 

Alternatively,  one  might  discard  the  ^-parameter  and  employ  a  modeling  scheme  similar  to  that  above, 
but  using  tower  measurements  along  with  the  ground  measurements  in  the  modeling  phase.  This 
procedure  might  be  most  powerful  if  a  and  8  are  the  entities  to  be  modeled,  given  i  fairly  restricted 
selection  of  possible  Earth  models.  A  question  in  this  case  is  whether  the  Earth  models  can  be  reasonably 
constrained  sufficiently  to  enable  reliable  estimates  of  a,  8,  and/or  a8. 

Another  possibility,  partially  explored  by  Christopher  Jekeli  and  the  author,  is  to  model  the  field 
above  the  surface  directly  from  the  boundary  values  on  the  surface  and/or  on  the  tower  without  reference 
to  an  Earth  model.  Whereas  avoiding  the  Earth  model  circumvents  questions  concerning  mass-distribution 
equivalence,  a  new  difficulty  arises  in  that  a  gauge  potential  enabling  workable  differential  equations  relies 
on  the  as-yet-unknown  Yukawa  field  in  its  boundary  conditions.  Thus,  the  price  of  more  efficient  and 
elegant  data  usage  appears  to  be  increased  difficulty  in  engineemg  such  a  scheme. 
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